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The neutrino-nucleon cross section is slightly larger than that for v — N. Therefore, v will escape 
more quickly from core collapse supernovae leaving the star v rich. A diffusion formalism is used to 
calculate the time evolution of the mu and tau lepton number densities. These quickly reach steady 
state equilibria. We estimate that a protoneutron star with a maximum temperature near 50 MeV 
will contain over 50 % more and v T than and v T . Supernovae may be the only known systems 
with large mu and or tau lepton numbers. 
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0^ ' Mu and tau neutrinos and antineutrinos are produced copiously in (core collapse) supernovae. These neutrinos 
with tens of MeV energies can not undergo charged current reactions. Therefore, it is assumed \i and r neutrinos and 
O !' antineutrinos have identical distributions. 

However, the cross section for z/-nucleon elastic scattering is somewhat larger then that for D — N scattering. This 
allows v to escape the star more easily leaving the supernova v rich. In this paper, we examine the effects of recoil 
and or weak magnetism corrections to the v — N cross sections. These corrections are of order E/M, where E is 
the neutrino energy and M the nucleon mass. They lift the degeneracy between neutrinos and antineutrinos. To our 
knowledge, all previous supernovae works (see for example [1]) assume equal v and v cross sections. 

We find a large excess of v over v in the star and a large mu and tau lepton number for supernovae. Possible 
\ implications of this for the neutrino signal, sensitivity to new physics, nucleosynthesis and other phenomena are 
0^ ' discussed at the end of this paper. A diffusion formalism is used for the neutrino transport. This is adequate for our 
purposes since most of the effect comes from the core of the protoneutron star well inside the neutrino sphere. 

We focus on mu and tau neutrinos. For these, the physics is very simple and clear. We explicitly discuss mu 
neutrinos. All of our results apply unchanged to v T . The effect is also present for electron neutrinos. Indeed, it will 

■ change the ratio of v e to v e and may produce significant consequences by changing the proton fraction. However, for 
u e there are also charged current reactions so the situation is more complicated. Therefore, we postpone a discussion 
of v e until the end of this paper. 

It is a simple matter to expand the v — N elastic cross section da/dfl to first order in E/M, see ref. [2] for example, 

,S : % = [<$(! + x) + c 2 a (3 - x)] [1 - 3^(1 - x)] ± 4c a (c„ + F 2 )^(l -x)}, (la) 

> ■ 

with x =cos# for scattering angle 6. Here the plus sign is for v and the minus sign for v. The vector (c v , F2) and 

■ axial (c a ) couplings are given in Table I. The transport cross section a — f dflda / dfl(l — x) is, 

2aiB \cl(l - 3^) + cl(5 2 A ± 8c a (c v + F a )£l . (16) 
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The weak anomalous moment F2 describes the a fJlU q 1 ' coupling of the Z to the nucleon. It could have significant 
strange quark contributions. Furthermore, the large value of F2 is important for this paper. A measurement of F2 
is underway using parity violating electron scattering [3]. For simplicity we neglect strange quark contributions to F2 
and c a in this paper. 

The difference in transport cross sections, to order E/M, is 

D _ Pv - ay = ScJfiv + F 2 ) E = E 

The coefficient 5 is 3.32 for neutrons and 2.71 for protons (using Table I). We assume that nucleon elastic scattering 
dominates the opacity, so D will also give the difference of mean free paths, D = (Ap — A„)/(A P + A„). 
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A simple diffusion equation for neutrinos is, 



^n(E) - Ul • A„Vn(£) = 0. (3) 



We subtract a similar equation for v to get, 



l[n(E) - n(E)] - • XV[n(E) - n(E)] = ~V • (AA)V[n(£) + n{E% (4) 



with A = (A p + X v )/2 = X Q E$/E 2 and 



Here the reference mean free path Aq is, 



AA = (Aj? — A„)/2 = XD. (5) 



A - 1 = 2^ (c , +5c 2 K5 (6) 

with E an arbitrary reference energy and p n the density of neutrons. For simplicity E/M terms are dropped in Eq. 
(6). Note that we assume pure neutron matter. A nonzero proton fraction should not change our results very much 
since 5 in Eq. (2) is similar for protons and neutrons. 

We integrate Eq.(4) over energy J d 3 E/(2ir) 3 n(E) — p, assume local thermodynamic equilibrium and work to lowest 
order in the chemical potential fi over the temperature T, 



to get, 



_-(TV) + V.A V^ y -V.A VT 2 . (7b) 



This equation describes the time evolution of the muon number density in the star. We also calculate the lepton 
number current, 

i.-*-^v & .-^n (8) 

The first term in Eq. (8) describes the conventional diffusion of lepton number while the second term comes from the 
diffusion of neutrino pairs which produce a lepton number current because of the difference in mean free paths AA. 
In steady state equilibrium J v — J„ = giving for the chemical potential 

a remarkably simple result. 

We now discuss numerical solutions of Eq. (7b) to bound the diffusion time for muon number in a supernova. To 
our knowledge there have been no previous estimates of this time. We find that muon number diffuses faster than the 
thermal energy. This may be because the heat capacity is proportional to the large baryon number. Furthermore, a 
low energy v with a long mean free path can effectively transport muon number [4]. Since diffusion is fast, we will 
assume steady state equilibrium and use Eq. (9) for some later results. 

First, consider a uniform star of density p = 5 x 10 14 g/cm 3 and (baryon) mass M tot = 1.5 Solar masses. This 
has a radius R of 11.25 km. A lower density surface will only speed up diffusion. Likewise we neglect nucleon Pauli 
blocking corrections to Ao in Eq. (6). These will also increase the mean free path. Perhaps the most extreme case is 
for the center of the star to be hot. Muon number must diffuse all the way to r=0. Therefore we consider a simple 
temperature distribution characteristic of the later stages of the protoneutron star cooling (perhaps 10 sec. after 
collapse [1]), 

T{r) = T (l - M(r)), (10) 
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with M(r) the enclosed mass divided by M tot . For a uniform density M = (r/R) 3 . We neglect the small time 
dependence of the temperature during the short simulation. 

We start from the initial condition /i = inside the star and adopt a somewhat arbitrary surface boundary condition 
that p is in equilibrium, given by Eq. (9), at the surface for all times. Our results are not very sensitive to this choice. 
Figure 1 shows /i as a function of time for a central temperature To of 35 MeV. Muon number diffuses into the center 
of the star so that eventually p(r = 0) rises to its equilibrium value given by Eq. (9). This rise happens quickly with p 
reaching half of its equilibrium value by t « 0.25 seconds^. In comparison it takes several seconds for the temperature 
distribution to change from surface peaked to one with a maximum at r — [5]. We conclude that lepton number 
diffusion is fast and expect p to be near its equilibrium value, Eq. (9). 

Equation (9) implies a density asymmetry, 

P^1 = ^»=^^±T = 498- (11) 
p + p 9C(3)T 54C(3)Af M' K ' 

This is large, 0.186 at T = 35 MeV. Equations (9) and (11) can be understood as follows. The density of v rises above 
that for v until the larger v density times a shorter mean free path balances pAp for antineutrinos. 

We integrate p — p over the star assuming Eq. (9) and the temperature profile in Eq. (10). This gives, 

Here To is a characteristic maximum temperature and the profile factor fi = 15 J r 2 dr / i? 3 [T(r)/To] 4 is one for Eq. 
(10) and is expected to be of order unity for other temperature profiles. We also calculate the total number of muon 
neutrinos, 

N + N= [ d Mp + P) = ^(^) 3 h, (13) 



with the profile factor f 2 = 12 L r 2 dr / R 3 [T(r)/T ] 3 . For the temperature distribution in Eq. (10), fi = 1. 
For example at To = 50 MeV we have, 

N - N = 9.40 x 10 53 , (14) 

N + N = 4.43 x 10 54 , (15) 

or N = 2.69 x 10 54 and N = 1.74 x 10 54 . This is a remarkably large asymmetry. The star contains 54 % more than 
Dfj,. One might expect E/M correction terms to be small. However, the coefficient 5 is large and the temperature, 50 
MeV, is high. We note that the asymmetry depends only linearly on To, see Eq. (11). Thus a decrease in To will not 
decrease the asymmetry greatly. 

Equation (14) is the muon lepton number of the supernova. Furthermore, it also gives the r number. We predict 
that supernovae are the only known systems with large p and or r number. This could impact new physics. For 
example, matter enhanced v oscillations could depend on the v density [6] . A weak long range force that couples only 
to /i or r number is difficult to observe. One may be able to use supernovae to set limits on such forces. 

We have found a large asymmetry for neutrinos recoiling against heavy nucleons. It is possible that the protoneutron 
star undergoes a transition from hadronic to a pure or mixed quark matter phase which could also contain strange 
quarks [7]. In a quark phase, one expects recoil corrections not of order E/M but E/m q where m q « M/3 is a 
constituent quark mass. This might produce even larger separations of v from v which could lead to an observable 
signature of quark matter. 

Perhaps the simplest model of quark matter is to assume the neutrinos scatter from nearly free quarks and that the 
structure of the quarks are simple without anomalous moments. Evaluating Eq. (2) for a down or strange quark using 
the couplings in table I gives, D w 1.01T/(M/3) = 3.03T/M. This is comparable to the neutron value D = 3.32E/M. 
The absence of the factor F2 reduces the quark value and approximately canceles the enhancement from the small 



1 Steady state equilibrium is reached even faster for a surface peaked temperature distribution. 
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mass. Therefore, we do not expect a strong sensitivity to a transition to quark matter. However, this discussion does 
emphasize that the separation between neutrinos and antineutrinos is occurring deep in the protoneutron star and 
may be sensitive to the properties of dense matter. 

We now discuss many-body corrections to Eq. (2) assuming a hadronic phase. The effects on the nucleons of 
relativistic kinematics, Fermi motion, Pauli blocking and nuclear mean ficlds[8,9] have been studied in ref. [10]. One 
might think that M in Eq. (2) will be replaced by the Dirac mass M* < M increasing the asymmetry. However, the 
anomalous moment term F?, is defined with the free nuclcon mass, see Eq. (28) of ref. [11]. Therefore its contributions 
do not increase. We find that D in Eq. (2) is almost unchanged in the medium. 

Burrows and Sawyer [12] argue that RPA correlations will greatly reduce cross sections in dense matter. However, 
fully relativistic RPA calculations using an interaction that is consistent with the equation of state give smaller 
corrections [10,13]. Furthermore, RPA corrections should effect v and v in about the same way. Therefore, we do not 
expect significant RPA changes to the cross section difference of Eq. (2). 

We now discuss the neutrino signal from a supernova. We divide the time into three periods. For a short initial 
time interval, the star is not in steady state equilibrium. During this period, which we estimate may only last for 
tens to hundreds of msec, the star radiates significantly more v than v. Next, the star will radiate equal numbers of 
v and v in steady state equilibrium for a relatively long period. Finally, as the star cools it must radiate away its 
lepton number. Therefore, there will be an ending period where more v are radiated than v. However, the fractional 
difference between v and D during this time is expected to be small because the star cools slowly. 

The number of neutrinos in the star at any one time, Eq. (15) is much less than the total number of muon neutrinos 
radiated (of order 10 57 ). Therefore, the first phase, with its excess P, can only last for a fraction of the total time. 
One should consider the possibility of observing the excess v during the initial phase. However, it may be difficult. 
Future work should calculate the neutrino signal in more detail. 

If the star undergoes a prompt collapse to a black hole, it will take the large lepton number of Eq. (14) with it. 
Assuming lepton number does not couple to a long range force, it will simply be lost down the black hole. However, 
the brief neutrino signal which precedes a prompt collapse could be significantly antineutrino rich. Depending on the 
time of the collapse, neutrinos from the more symmetric second and third periods may never escape the star. 

Finally, we consider electron neutrinos. The neutral current effect that we have calculated for and v T will also 
apply to v e . In addition, charged current reactions have similar weak magnetism corrections. These terms will increase 
the v e cross section with respect to v e . Thus, the number of v e compared to D e will increase in the star. This could 
have a significant effect on the dynamics of the explosion and or nucleosynthesis by increasing the proton fraction. 
We will discuss this further in a later work [5] . 

The separation between neutrinos and antineutrinos involves parity violation. If the weak interactions conserved 
parity, the v and v cross sections would be equal (to lowest order in Gf)- Thus, the large mu and tau lepton numbers 
of a supernova can be considered a macroscopic manifestation of parity violation. 

In this paper we have calculated the effects of recoil or weak magnetism corrections to v — N cross sections. These 
are of order the neutrino energy E over the nucleon mass M, E/M and make the v — N cross section larger than that 
for v — N . A diffusion formalism was used to follow the time evolution of lepton number. It quickly reaches a steady 
state equilibrium where a larger v density compensates for the longer v mean free path. For a maximum temperature 
near 50 MeV, we estimate that protoneutron stars contain over 50 % more and v T than and v T . Core collapse 
supernovae may be the only known systems with large mu and or tau lepton numbers. 
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TABLE I. Vector and axial couplings for neutral current scattering from neutrons, protons and up and down quarks. We 
assume sin 2 (9=0.23 and g a — 1.26. 
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FIG. 1. Muon neutrino chemical potential vs. radius. The temperature distribution is assumed to be given by Eq. 
(10) with a maximum central temperature of To = 35 MeV. The curves are for times of 0.02 (dotted) to 1 second 
(dot-dashed). The solid curves are for intermediate times of (bottom to top) 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6 seconds. 
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